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ABSTRACT 

We present the first > 3.5cr (conservative) or > 5.8cr (sum of lines significance) 
detection of two Warm-Hot Intergalactic Medium (WHIM) filaments at 2; > 0, 
which we find along the line of sight to the blazar Mkn 421. These systems 
are detected through highly ionized resonant metal absorption in high quality 
Chandra-AClS and -HRC Low Energy Transmission Grating (LETG) spectra of 
Mkn 421, obtained following our two Target of Opportunity requests during two 
outburst phases (-Fo.5-2fcey = 40 and 60 mCrab; 1 Crab = 2xl0~® erg s~^ cm~^). 

Columns of He-like oxygen and H-like nitrogen can be detected in the coadded 
LETG spectrum of Mkn 421 down to a sensitivity of Nqvii > 8 x 10^"^ cm~^ and 
^Nvii > 10^^ cm^^ respectively, at a significance > 3a. The two intervening 
WHIM systems that we detect, have OVII and NVII columns of Nov// = (1.0 ± 
0.3) X 10^5 cm-2 Nnvii = (0.8 ±0.4) x 10^^ cm'^, and Nqvii = (0.7 ±0.3) x 10^^ 
cm~^, l^Nvii = (1-4 ± 0.5) X 10^^ cm~^ respectively. We identify the closest of 
these two systems with an intervening WHIM filament at cz = 3300 ± 300 km 
s~^. The second system, instead, at cz = 8090 ±300 km s~\ is identified with an 
intervening WHIM filament located ~ 13 Mpc from the Blazar. The filament at 
cz = 3300 ±300 lies ^5 Mpc from a known HI Lya system at cz = (3046 ±12) km 
s~^ (Shull et al., 1996) whose 3a maximal HI kinetic temperature, as derived from 
the observed line FWHM, is T < 1.2 x 10^ K. This temperature is inconsistent 
with the temperature measured for the X-ray filament, so if the systems are 
related a multiphase WHIM is required. 

Combining UV and FUV upper limits on the HI Lya and the OVl2s^2p tran- 
sitions, with our measurements in the X-rays, we show that, for both filaments, 
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equilibrium coUisional ionization (with residual photoionization by both the dif- 
fuse UV and X-ray background, and the beamed emission of Mkn 421 along our 
line of sight) provides acceptable solutions. These solutions define ranges of tem- 
peratures, metallicity ratios and equivalent H column densities, which are in good 
agreement with the predictions of hydrodynamical simulations for the formation 
of large scale structures in the Universe. 

Prom the detected number of WHIM filaments along this line of sight we 
can estimate the number of OVII filaments per unit redshift with columns 
larger than 7 x 10^^ cm'^ dMovii/dz{Novn > 7 x 10^^) = 67l43, consistent, 
within the large la errors, with the hydrodynamical simulation predictions of 
dMovii / dz{Novn > 7 x 10"^^) = 30. Finally, we measure a cosmological mass 
density of X-ray WHIM filaments fi^ = 0.0271°;°?^ x lOt^/^l-i, consistent with 
both model predictions and the estimated number of 'missing' baryons at low 
redshift. 

Subject headings: 

1. Introduction 

Despite dramatic recent progress in cosmology, very little is yet known about the location 
of the normal baryonic matter in the local Universe. The extraordinary WMAP ( Wilkinson 
Microwave Anisotropy Probe) measurements of the Cosmic Microwave Background (CMB) 
anisotropies favor a A-CDM scenario in which the baryon density in the Universe amounts 
to Qf, — (4.6 ± 0.2)hjQ % of the total matter-energy density, with a baryon-to-matter ratio 
of 0.17 ±0.01 (Bennett et al., 2003, Spergel et al., 2003). This number agrees well with pre- 
dictions by the standard 'big bang nucleosynthesis' when combined with light element ratios 
{rib = (4.4 ± OA)hjQ %: Kirkman et al., 2003), and also with the actual number of baryons 
detected at redshifts larger than 2 in the 'trees' of the Lya Forest: > 3.5hjQ % (Ranch, 
1998; Weinberg et al., 1997). This concordance cosmology represents a major advance. How- 
ever, the number of baryons actually detected in the virialized Universe (i.e. stars, neutral 
H and He and molecular H associated with galaxies. X-ray emitting gas in clusters and big 
groups with M> 4 x 10^^ M©: Fukugita, 2003) and the non yet virialized Universe (i.e. pho- 
toionized warm/cool gas in the local Lyman-a forest filaments and shock- heated warm-hot 
gas in the OVI filaments: ShuU, 2003; Stocke, ShuU & Penton, 2004; Danforth & Shull, 
2005), is less than half of the expected concordance value, Qi,{detected) = (2.1 ± 0.3) % (for 
/i7o = 1). The number of 'missing baryons' at ^$2 is therefore flb{missing) = (2.5 ± 0.4) %. 

If we cannot find half of the ordinary matter in our local Universe it does not bode well 
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for understanding dark matter or dark energy. 

Hydrodynamical simulations for the formation of structures in the Universe, offer a 
possible solution to this puzzle of the "missing baryons" . They predict that in the present 
epoch {z^l — 2), 30-40 % of the ordinary, baryonic, matter in the Universe lies in a tenuous 
medium at overdensities 5 ~ 5 — 200 ^ , relative to the mean baryon density in the Universe 
(~ 10^ below the density of the interstellar medium in our Galaxy). This medium is in a 
"warm-hot" phase and forms the so-called 'Warm Hot Intergalactic Medium' or "WHIM" . 
The WHIM was shock- heated to temperatures of 10^ — 10^ K during the continuous process of 
structure formation (e.g. Hellsten et al., 1998; Cen & Ostriker, 1999; Dave et al., 2001; Fang, 
Bryan & Canizares, , 2002). At these temperatures C, N, O and Ne (the most abundant 
metals in gas with solar or solar- like composition), are mainly distributed between their He- 
like and H-like species. The strongest bound-bound transitions from these ions fall in the 
soft X-ray band (e.g X{OVIlKa) = 21.602 A, and X{OVIIlKa = 18.97 A), and so a 'forest' 
of weak (i.e. Nx< < 10^^ cm~^, or Wxi < 18 — 30 mA depending on the particular ion ^ 
) metal absorption lines is expected to be imprinted in the X-ray spectra of all background 
sources at redshifts between zero and ~ 1 — 2 , by the WHIM filaments (e.g. Perna & 
Loeb, 1998; Hellsten et al., 1998; Fang, Bryan & Canizares, 2002). These lines are the direct 
analogous of the 'Lyman-a forest' due to cool HI and detected copiously at higher redshift 
in the optical spectra of quasars. 

Ubiquitous high- velocity OVI absorption (Sembach et al., 2003) and higher- ionization, 
CV-VI, NVI-VII, OVII-VIII and NelX-X absorption (Nicastro et al., 2002, Fang et al., 2003, 
Rasmussen et al., 2003, Cagnoni et al., 2004, Williams et al., 2005) at 2; ~ have been 
detected in the far-ultraviolet (FUV) and X-ray spectra of extragalactic sources. The exact 
location of this ubiquitous absorber, whether in an extended Galactic corona (Sembach et 
al., 2003, Gollins et al., 2004) or in a local multiphase WHIM filament containing our own 
Local Group of Galaxies (Nicastro et al., 2002, Nicastro et al., 2003, Williams et al., 2005) 
is still highly controversial. Either case the success in detecting such a diffuse medium at 
2; = is due to our special embedded location and so on the conspicuous portion of absorber 
crossed by any single extragalactic hue of sight. 

Detecting the WHIM outside the Local Group, instead, means observing along a ran- 
dom, unprivileged line of sight. Due to the intrinsic steepness of the 'LogA/'-Log(N^)' of 
the WHIM metal absorption hues, at relatively large columns (e.g. Hellsten et al., 1998, 



= {rib/ < n-b >) is the baryon gas overdensity in units of average density in the Universe, < Ub >= 
2 X 10-^(1 -I- zfi^bh"^ /0.02) cm-3. 

^here, and throughout the paper, is the generic He- or H-Uke ion of C, N, O or Ne 
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Fang, Bryan & Canizares, 2002), only a single OVII Ka system with Nov// > 10^^ cm~^ is 
expected along a random line of sight, up to 2; = 0.3, with a probability of ~ 60 %. This 
number rises to 8 for systems with Nqvii x 10^^ cm~^. So detecting weaker lines is expected 
to be more productive than probing higher redshifts. 

Consistently with expectations, no secure X-ray detection of 2; > intervening WHIM 

absorption systems has been claimed so far. The deepest observation of a quasar at 2; ~ 0.3 
is the 500 ks Chandra-LETG observation of H 1821+643. This observation is sensitive only 
to N^i t7 X 10^"'' cm^^ (at > 3cr level), and only 2-3 low-significance (~ 1 — 2a) detections of 
such systems have been claimed (Mathur, Weinberg, & Chen, 2003). Ravasio et al. (2005) 
recently published a moderate signal to noise Newton-XMM Reflection Grating Spectrometer 
(RGS) observation of the neraby {z — 0.03) of Mkn 421. This spectrum is sensitive to OVII 
columns of ~ 3 — 4 x 10^^ cm~^, at 3a and no intervening systems could be identified down 
to this OVII column sensity threshold. Finally, another 2; > OVIII WHIM candidate, 
at z — 0.0554 was claimed by Fang et al. (2002), in a Chandra Advanced CCD Imaging 
Spectrometer (ACIS) LETG spectrum of the blazar PKS 2155-304. However neither a deeper 
C/ian(/ra High Resolution Camera (HRC) LETG nor a Newton-XMM observation of the same 
object could confirm this result (Nicastro et al., 2002, Rassmussen et al., 2003, Cagnoni et 
al., 2004). Clearly higher signal-to- noise spectra are urgently needed. 

A very large number of counts per resolution element (CPREs) is needed in current 
Chandra Low Energy Transmission Grating (LETG) or Newton-XMM Refiection Grating 
Spectrometer (RGS) spectra, both of which have a resolution of R~ 400 at ~ 20 A, to 
detect these lines. In an ideal detector about 30 counts, at 21.6 A, are needed in a AA = 50 
mA continuum bin (i.e. 1 resolution element) for a 3a detection of a 29 mA absorption line. 
In practice, at such a low number of CPREs and for moderately noisy detectors (like the 
Chandra HRCS-LETG), at least 100 CPREs are needed. Such a line probes OVII columns 
of about 10^^ cm~^, the very high-column tail of the WHIM distribution. To probe the 
more common WHIM columns of Nov 11 ~ 10^^ cm~^, about 2500 CPREs are needed for 
a 3a detection. The brightest Scyfcrt galaxies and blazars have typical soft X-ray (0.5-2 
keV) fiuxes of the order of ~ 4 x 10^^^ erg s^^ cm~^ (2 mCrab), which would take a 2.5 
Ms Chandra-LETG exposure to reach this S/N. Furthermore all these objects are nearby 
{z^0.05) and so sample only path-lengths of a few hundred Mpc. 

A possible solution is to observe background sources in untypically bright states. This 
strategy has been previously and successfully exploited with low resolution spectrometers 
(e.g. 5e^^oSAX) to study Gamma Ray Burst X-ray afterglow, but also compact Galactic 
sources and non-beamed AGNs (e.g. Puccetti et al., 2005, in preparation). We have extended 
this strategy to beamed, pure continuum, AGNs in order to obtain exceptionally good quality 
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spectra of the intervening Intergalactic Medium (IGM). In this paper we present the first 
results of this highly successful observational strategy. We report the first high significance 
(cr = 3.5 and a = 4.8) detections of two high ionization absorbers at cosmological distances, 
which we identify with WHIM filaments, in a high signal-to-noise Chandra-LETG spectrum 
of the blazar Mkn 421. A companion, complementary paper, which reports on the first 
estimate of the cosmological mass density of OVII WHIM filaments derived from these 
detections, has been recently published (Nicastro et al., 2005). 

Throughout the paper we use LETG spectra grouped at resolutions of AA = 12.5 mA 
(for fitting purposes) and 25 mA (for plotting purposes), i.e. respectively 4 and 2 times 
better than the intrinsic LETG resolution of 50 mA. We adopt Hq = 72 km s^^ Mpc^^ 
(Preedman et al., 2001, Bennett et al., 2003). Errors are quoted at la for 1 interesting 
parameters, unless otherwise stated. 

2. The Chandra Data of Mkn 421 

Mkn 421 is a blazar at z = 0.03 (Dc Vaucouleurs et al., 1991), which undergoes frequent 
strong outbursts from its quiescent state of ~ 1 mCrab (2 x 10~^^ erg s^^ cm^'^) to 20- 
50 mCrab (Figure 1). When in outburst Mkn 421 is one of the few objects detected at 
TeV energies (e.g. Okumura et al., 2002). Mkn 421 has been observed by Chandra at four 
different epochs (Figure 1, Table 1). The first short observation was performed with the ACIS 
High Energy Transmission Grating (HETG) configuration, at the beginning of the satellite 
calibration phase, on 1999 November 5. This observation lasted 25.8 ks, and the target was 
caught at a historical minimum in its 0.5-2 keV fiux, of -Fo.5=2feeF = 10"^^ erg s^^ cm~^ = 0.5 
mCrab, or 0.75 mCrab extrapolating to the 1.5-12 keV band of the Rossi-XTE All Sky 
Monitor (ASM: e.g. Levine et al., 1996: Figure 1). The number of photons per resolution 

Table 1: Table of X-ray Observations 



Date 


MJD 


Det/Grating 


Exposure 


(0.5-2) keV 


Total Dispersed 


CPREs 










Flux 


Counts 


at 21 A 








ks 


mCrab 






11/05/99 


51487 


ACIS/HETG 


25.8 


0.5 


9635° 


2« 


05/29/00 


51693 


ACIS/HETG 


19.6 


9.5 


125301" 


17" 


05/29/00 


51693 


HRC/LETG 


19.7 


13.5 


209586'' 




10/26/02 


52573 


ACIS/LETG 


91.5 


60.0 


4213181'= 


2676" 


07/01/03 


52821 


HRC/LETG 


99.2 


40.0 


2967358^ 


2512'' 



" 1st order MEG. All orders LETG. 1st order LETG. 
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elements (AA = 20 mA) in the Ist-order HETG spectrum of Mkn 421 was therefore far below 
the threshold needed for detection of WHIM absorption lines, and so useless for our purposes 
(Table 1). We no longer consider this observation here. Six months later (2000, May 29) we 
requested and obtained two short consecutive Chandra Director Discretionary Time (DDT) 
observations of Mkn 421, following 2 weeks of intense activity of the target as recorded by 
the ASM (Figure 1). The two observations were performed ~ 1 week after the trigger with 
the Chandra ACIS-HETG and HRC-LETG, and lasted ~ 20 ks each, recording 0.5-2 keV 
fluxes some 20 times higher than before, 9.5 and 13.5 mCrab respectively (Figure 1, Table 
1). Despite the large fluxes, only 177 and 17 CPREs were detected at A = 21.0 A in the 
LETG (all-orders) and HETG (Ist-order) spectra, respectively (Table 1). In the following we 
will use only the HRC-LETG spectrum and only in combination with the subsequent, much 
better quality, LETG spectra of the source. Finally, two further 100 ks Chandra ACIS-LETG 
and HRC-LETG observations of Mkn 421 were performed in 2002, October, 26-27 (Figure 1 
and 2, Table 1), and 2003, July 1-2 (Figure 1, Table 1), under our approved Chandra-KOA 
program to observe blazars in outburst phases. These two observations were performed ~ 24 
and ~ 18 hours respectively after the trigger and caught the source at a historical maximum, 
~ 100 times brighter than in the first observation, 60 and 40 mGrab in the 0.5-2 keV band 
(Figure 1, Table 1), and allowed us to collect ~ 2600 GPREs each at 21 A (Table 1). 



3. The Combined LETG Spectrum of Mkn 421 

To increase the signal to noise, we combined the first order ACIS-LETG spectrum (and 
responses) of Mkn 421, with the two non-order-sorted HRC-LETG spectra (and responses) of 
the target in the overlapping 10-60 A wavelength range. The final co-added LETG spectrum 
contains more than 7 miUion counts, and has ~ 5300 CPREs in the continuum at 21 A. This 
is the best S/N of any high resolution spectrum taken with Chandra, including Galactic 
X-ray binaries. This S/N was enough to detect OVII columns of Novate x 10^*^ cm~^ at a 
significance level > 3cr. 



3.1. Spectral Responses and Continuum Fitting 

ACIS-LETG negative and positive 1st order effective areas (ARFs) were built by using 
the Jullgarf script distributed with Ciao (vs. 3.0.2), with the ACIS-Quantum- Efficiency 
(QE) file made available by the ACIS calibration group (CALDB 2.2.6) to correct for the 
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Fig. 1. — 1.5-12 keV Rossi-XTE ASM lightcurve of Mkn 421, from 1999, July 29 to 2004, April 14. 
Superimposed starred points are 0.5-2 keV Chandra flux measurements extrapolated to the 1.5-12 keV ASM 
band by using the best fitting Chandra continua. 



- 8 - 





Fig. 2. — Portion of the 2-D ACIS-LETG dispersed spectrum of Mkn 421 during our first Chandra TOO 
on 2002, October, 26-27. The dispersed spectrum is seen in the two very bright hues which form diagonal 
from the upper left to the bottom right, the positive and negative orders, and this is the highest signal-to- 
noise spectrum taken so far with the LETG . The "fan of rays" around both the Oth-order image and the 
dispersed spectrum are an instrumental artifact due to diffraction of the X-ray photons off the instrument 
support. The lack of photons (i.e. the "hole") in the very center of the Oth-order source image is due 
to an exceptionally high degree of "pile-up" in the CCD detector. Finally the luminous strip centered on 
the Oth-order image perpendicular to the dispersed spectrum contains Oth-order photons deposited during 
read-out along the CCD read-out axis. This is because the source photon arrival rate is much faster than 
the 0.8 s CCD read-out rate. 
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ACIS contamination problem ^ (e.g. Marshall, et al., 2003). Negative and positive ACIS- 
LETG 1st order ARFs were then checked for consistency and combined using the Ciao tool 
add_grating_orders, and finally convolved with the standard 1st order ACIS-LETG Redistri- 
bution Matrix (RMF, CALDB 2.2.6). For the HRC-LETG responses, instead, we used the 
standard convolution products of RMFs and ARFs files, for orders 1 through 6 ^ . Finally, 
the first-order ACIS-LETG and order 1 through 6 HRC-LETG responses were combined 
together (weighting by count rate) to form a ACIS-LETG-I-HRC-LETG response matrix. 

We used the fitting package Sherpa (Freeman et al., 2001), in Ciao, to model the broad 
band 10-60 A ^ average continuum of Mkn 421 during the three LETG observations with a 
single power law absorbed by neutral Interstellar matter with solar composition (Morrison 
& McCammon, 1983) along the line of sight. The best fitting hydrogen column of Nh = 
(1.41 ± 0.01) X lO^'' cm~^ is in excellent agreement with the Galactic value along this line of 
sight (iVf"' = 1.4 X 10^0 cm-^: Elvis, Lockman & Wilkes, 1989; Dickey & Lockman, 1990), 
and the best fitting power law has a slope of F = 1.802±0.003. This model, however, does not 
yield an acceptable fit to the data (x^(d.o./.) = 3.8(1635)). While the large scale continuum 
shape is well accounted for, there are instead highly significant small scale systematic positive 
and negative residuals across the entire wavelength band (Figure 3). Hence attempts to 
reduce the differences between the model and the data by the standard method of introducing 
an additional continuum component, as a second power law, a broken power law, or a black- 
body, do not yield satisfactory results. This is due to the localized nature of the strong 
deviations, even the largest of which (at ~ 43 A) would fit into a single resolution element 
of earlier CCD spectra. 



3.1.1. Isolating Instrumental Artifacts in the Spectrum 

The strongest negative and positive features in the ratio between data and best fitting 
model, fall close to the neutral C and O K-edges respectively, at ~ 23 A and ~ 43 A (Fig. 3), 
as well as some of the ACIS-LETG bad columns associated to node-boundaries ^ or chip-gaps 

■^http: / / asc.harvard.edu/ciao/why/acisqcdcg.html 

''http: / /asc. harvard. cdu/cal/Lctg/Hrc_QE/ca_indcx.htnil:^rsps 

^We also excluded the wavelength range 48-57 A, where the two HRC plate gaps are present, and non 
modeled in the effective area 

^ACIS Node-boundaries (3 per chip) are normally associated with bad columns. These columns, however, 
are not completely blind in wavelength scale, due to the satellite dithering, and appear as smooth and broad 
1 A) notch-like features in the aspect-reconstructed dispersed spectrum. 
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( particularly, the ones at ~ 20 A, and ~ 13 A, Fig. 3). All these features are most likely 
due to residual calibration uncertainties both in the ACIS- and HRC-LETG. In particular 
residuals from the CI and 01 K-edge wavelengths shortward are likely to be mostly due to 
uncertainties in the model description of the instrumental absorption caused by a layer of 
contaminant material that has gradually condensed on the ACIS window, and that affects all 
observations taken with the ACIS detector (Marshall et al., 2003) ^ . Clearly the depth of the 
CI edge is severely underestimated by the model (amplitude of the data/model ratio ^ — 30 
%, Fig. 3), while the 01 edge depth is overestimated (amplitude < +10 %, Fig. 3). The exact 
shape of the recovery of the continuum level shortward of these two edges is important in 
setting the proper level of the local continuum, and so to properly evaluating the significance 
and amplitude of weak narrow absorption features in these spectral regions. As an example. 
Figure 4a shows the 20-25 A region of the LETG spectrum of Mkn 421, with the best fitting 
powerlaw plus neutral absorption continuum superimposed. The fit is unacceptable, with a 
reduced-x^ of Xr{d-o.f.) = 6.98(197). Clearly the 01 K-edge is overestimated by the model, 
and so is all of the continuum between ~ 21.7 and 25 A, leaving strong systematic residuals 
in the data/model ratio. The strong narrow absorption features in the spectrum are OVII 
and 01 absorption a,t z — (Fig. 4a). We also note here the scatter of the data in the 
narrow portion of the spectrum immediately longward of the OVII Ka line at z — and up 
to ~ 22.3 A (where the OVII Ka transitions at the rcdshift of Mkn 421 would fall), which 
may suggest the presence of additional narrow line-like features. Figure 4b shows the highly 
non-gaussian (60 bins outside the (-2,+2)o" interval, versus the expected 7.5 in the gaussian 
hypothesis), positive-skewed (85 bins fall in the (0,+3)o' interval, versus the 55 bins in the 
corresponding negative interval) and apparently bimodal (reflecting the systematics in the 
data/model ratio) distribution of the residuals in a in relative number of bins per half-o", 
after ignoring narrow intervals of data at the wavelengths of the strong z — OVII, atomic 
01 and molecular 01 absorption fines. 

We thus decided to re-fit the broad-band data leaving the metal abundances of the 
photoelectric absorption model used to model the ISM absorption, free to vary. This fit gave 
a best fitting H absorbing column of (1.44 ±0.01) x 10^° cm"^ (again in agreement with the 
Galactic ISM column observed in neutral H, Elvis, Lockman & Wilkes, 1989), with C and 
Ne over-abundance of 3 and 4 times solar respectively, and N and O virtually absent. We 
stress that these values should not be thought of as true differences in the ISM metallicity, 
compared to solar, but rather as due to calibration imcertainties. The best fitting power 
law slope now becomes ~ 0.1 steeper, F = 1.911 ± 0.003. The continuum model yields 
much better descriptions of the CI and 01 K-edges (Fig. 5a,b), but it slightly amplifies the 



''http: / / cxc.harvard.edu/cal/Acis/CaLprods/qeDeg/ 
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Fig. 3. — Ratio between the combined LETG spectrum of MKN 421 and its best fitting simple power 
law absorbed by neutral ISM with metal abundances fixed to solar value. NB = 'Node Boundary'; FB = 
'Front-Illuminated to Back-Illuminated'; BF — 'Back Illuminated to Front-Illuminated'. 




Fig. 4. — (a) 20-25 A region of the LETG spectrum of Mkn 421 and best fitting powerlaw plus neutral 
absorption continuum (solid curve), (b) Distribution of the fraction of residuals per half-cr with a. 
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discrepancies between data and model at the wavelengths of several ACTS node-boundaries 
or chip-gaps, particularly around 13 and 20 A (Fig. 6a,b). Moreover locally the continuum 
level shortward of the CI and 01 K-edgcs continues to be somewhat over /under-estimated 
(Fig. 5a,b). This may at least partly reflect the fact that the photoelectric model we use to 
describe the ISM absorption uses photoelectric cross-sections that do not properly describe 
the complicated absorption structure longward of the instrumental edge wavelengths * . To 
allow for these deviations we added ten broad gaussians (6 in emission and 4 in absorption). 
These have best fitting FWHM of 0.5-7 A, except at the exact wavelength of the CI K-edge 
where a relatively narrow gaussian, with FWHM of 0.2 A fits well. With this adjustment 
we recover a good fit to the continuum level in three important regions of the spectrum 
longward of the CI edge (between ~ 40 — 43 A, ~ 26 — 30 A and ~ 31.5 — 35.5 A), and at 
the wavelengths of two ACIS node-boundaries (~ 12 — 14 A, and ~ 19 — 21 A). This model 
yielded a Xri^'^f) ~ 1-81(1598), and left unchanged both F and Nh- Fig. 7a shows the 
residuals, in a, between the data and our final best fitting continuum model. The residuals 
are now relatively fiat over the entire 10-48 A wavelength range (Fig. 7a). There are still 
relatively broad excesses at the 01 and CI K-edges wavelengths. There are also several 
narrow negative spikes, which we show in the next section, all due to intervening narrow 
resonant absorption metal lines. Fig. 7b show the normalized distribution of residuals to 
our best fitting continuum model, compared with a gaussians distribution. The positive tail 
of the distribution is gaussian. The negative-only < —3a tail in the distribution is due to 
the presence of several real and high significance absorption lines in the data. 



3.2. Intervening Resonant Absorption Lines 

Our statistical analysis (Figure 7b) shows that an excess, compared to a gaussian dis- 
tribution, of 28 > 3a negative residuals, with 16 of these being > 4(7, are present in the 
LETG spectrum of Mkn 421, relative to the best fitting continuum model. Inspection of 
relatively narrow regions (AA = 2 — 5 A) of the spectrum (Fig. 8a, f), confirms these findings 
and shows that several groups of 2-4 contiguous 25 mA negative bins deviate from the best 
fitting continuum with significance > 1.5a per bin (i.e. > 3a per group). It is striking that 
these narrow absorption line-like features concentrate in critical spectral regions, always at, 
or longward of, the rest-frame positions of the main resonant transitions from highly ionized 
or neutral C, N, O and Ne (Fig. 8). We therefore tentatively identified these features as metal 
absorption lines from either the interstellar medium (ISM) of our Galaxy, the Local group 
WHIM Filament (Nicastro et al., 2002, Nicastro et al., 2003, WiUiams et al., 2005), interven- 



'http://space.mit.edu/CXC/calib/letg_acis/ck_cal.html 
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Fig. 5. — LETG spectrum of Mkn 421 and best fitting powerlaw plus neutral absorption with varying 
metallicities continuum (solid curves) around (a) the CI K-edge and (b) the 01 K-edge. 
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Fig. 6. — LETG spectrum of Mkn 421 and best fitting powerlaw plus neutral absorption with varying 
metallicities continuum (solid curves) around (a) 19-20 A and (b) 12.8-13.8 A ACIS node-boundaries. 
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Fig. 7. — (a) Residuals in a of the LETG spectrum of MKN 421 to its best fitting continuum model: a 
power law absorbed by neutral ISM with metal abundances free to vary, plus ten broad gaussians used to 
recover the continuum level in spectral regions still affected by non-optimally modeled instrumental features. 
Rest Frame wavelengths of the Ka transitions from He-like and H-like C, N, O and from neutral O, are 
labeled, (b) Normalized distribution of residuals to the best fitting continuum model described above (solid 
histogram), compared with the gaussian distribution (dashed histogram). 




Fig. 8. — Six portions of the LETG spectrum of Mkn 421 along with it best fitting continuum plus narrow 
absorption model (solid line), centered around the rest-wavelengths of the (a) NeXI — Xxa: (b) OYlllxa, 
(c) OVllfCa, (d) NVIIxq, (e) NVlxa and (f) CVI^q transitions. For each labeled ion, the three vertical 
lines from left to right are the rest frame wavelength (thin line) and the expected wavelengths at z = 0.011 
(medium thickness line) and z = 0.027 (thick line). 
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ing intergalactic absorbers, or material outflowing from the nucleus of Mkn 421 along our line 
of sight. Wc fitted these features with negative gaussians, leaving all parameters (FWHM, 
centroid and amplitude) free to vary. Our best fitting model thus contains a total of 24 un- 
resolved (i.e. intrinsic FWHM < 50 mA) absorption lines, and has a Xri'^'^f) ~ 1-25(1529) 
(giving an F value of 14.2, when compared with the best fitting continuum model with no 
lines, corresponding to a probability of exceeding F, P(> F) « 0.001). Figure 9 shows 
the normalized distribution of residuals to our best fitting model, including absorption lines, 
over the 10-48 A range, compared with a gaussian distribution. The negative tail present 
in the distribution of Fig. 7b, is no longer visible, and the overall distribution is now close 
to gaussian. We are therefore confident that the measured significance of our absorption 
lines, as expressed in number of standard deviations from a parent gaussian distribution, 
represents a good estimate of the significance of the lines. 

Figure 8 (a to f) shows 21 of the 24 absorption fines fitted ^ in 6 different portions of 
the LETG spectrum of Mkn 421, centered around the rest- wavelengths of the NeXI-X/<:a, 
OVIIIxa, OVllKa: N Vilify, NVI^a and CVIxa transitions. 

3.2.1. Line Identification 

We identify 14 of the 24 absorption lines in the LETG spectrum of Mkn 421 as being 
produced by 2; = ± 0.001 {cz — ±300 km s~^) systems of neutral (molecular and atomic 
01) and ionized metals. Physically the system corresponds to two distinct components: the 
ISM of the Milky Way (i.e. Savage et al., 2000, Savage et al., 2003, Savage et al, 2005, 
ApJ) and the Local Group WHIM (e.g. Nicastro et al., 2002; Nicastro, 2003, and references 
therein). One more line, at A = (22.32 ± 0.02) A (Figure 8c) is tentatively identified here 
with the combination of the z = OV Ka inner shell transition {/^rest-frame = 22.37 A 
Schmidt et al., 2004) and an OVII Ka transition at z = 0.033 ±0.001 (~ (900 ±300) km s"^ 
bluer than the estimated redshift of the blazar's host galaxy: De Vacouleurs, 1991). These 
local systems are the subject of a forthcoming paper (Williams et al., 2005), and will not be 
discussed further here. 

The remaining 10 lines could not be identified with any other z — Q transition from either 
neutral or ionized metals. Two fines (A = 13.80 A and A = 25.04 A, see Table 2) fall close to 
the rest frame wavelengths of two outer shell transitions from FeXVII, at A = 13.823 A, and 
ArXII at A = 25.04 A (Verner, Barthel & Tytler, 1994). However, these two transitions have 
oscillator strengths [f {FeXVII) = 0.331 and f {ArXII) = 0.191: Verner, Barthel & Tytler, 



^Not shown are the neutral O from the ISM, and the CV Ka at ^; = (WilUams et al., 2005). 
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Fig. 9. — Normalized distribution of residuals to our best fitting model, including 24 narrow (FWHM< 50 
mA) absorption gaussians. The dashed histogram is the gaussian distribution. 
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1994] that are respectively, 8.9 and 7.6 times weaker than the oscillator strengths of two 
other transitions from the same ions: the FeXVII at A = 15.01 A (/ = 2.95) and the ArXII 
doublet at < A >= 31.37 A (/i + /2 = 1.45). These are two regions of the spectrum not 
directly contaminated by known instrumental features and with effective areas comparable, 
or slightly lower, than the effective areas at the wavelengths of the detected features. So both 
the FeXVII;^i5 01 and the ArXlI;^3i 37 lines should be detected in this spectrum at a much 
larger significance than that of the putative z — hnes at 13.80 and 25.04 A respectively. 
No hne is instead detected at these wavelengths so ruhng out the z — FeXVII and ArXII 
identifications for these two features. Finally, although the vast majority of the inner shell. 
X-ray, transitions from C, N O and Ne are poorly known, both in wavelength and strength 
(ATOMDB, v.1.3.0. Smith et al., 2001; Behar & Netzer, 2002), basic atomic physics requires 
that the strongest of these transitions (Ka) must all fall at longer wavelengths than the Ka 
He-like transitions from the same elements. 

We detect 5 lines in the OVII Ka z = — 0.033 range (Figure 8c). From the left, 
the first line is unambiguously identified with the OVII Ka at 2; = 0. The next line to the 
right, at A = 21.85 A, cannot be due to inner-shell absorption by the lower-ionization ion of 
oxygen, OVl, since the strongest transition from this ion (ls^2s Is2s2p) falls at A = 22.03 
A (HULLAC value: Bar-Shalom et al., 2001; Liedhal, D.A., private communication), that is 
redward of the line at A = 21.85 A (moreover the OVI Ka line at 2; = is indeed detected in 
our spectrum, as labeled in Figure 8c; see also WiUiams et al., 2005). Similarly, the next hne 
to the right, at A = 22.20 A, cannot be due to OV absorption, since the strongest inner-shell 
transition from this ion falls at A = 22.37 A (Schimdt et al., 2004), that is redward of the 
line at A = 22.20 A. A weak OV Ka inner-shell transition may indeed be detected in our 
spectrum, in the red wing of the line at A = 22.32 A(Figure 8c). Finally, three strong (i.e. 
oscillator strength > 0.07) inner-shell transitions from the next lower ionization ion, OIV, 
have rest frame wavelengths A = 22.73 — 22.75 A(unresolved in the LETG: HULLAC values, 
Liedhal, D.A., private communication). These wavelengths are longward of the OVII Ka 
wavelength at the redshift of Mkn 421. We conclude that the two weak lines at A = 21.85 
and A = 22.20 A must be due to highly ionized intervening O absorption at 2; > 0. 

We repeated this analysis for the other three portions of the spectrum of Mkn 421 
that maybe in principle affected by inner-shell C, N and Ne absorption, redward of the 
Ka transitions from the He-hke ions of these elements. We concluded that only two of 
the remaining 8 features may in principle be due to inner shell transitions from moderately 
ionized Ne or N: the line at A = 13.80 A (Figure 8a) and the hne at A = 29.54 A (Figure 8e; 
see Table 2). The line at A = 13.80 A lies close to the rest frame wavelength of the forbidden 
line of the NelX triplet, analogous to the OVI Ka inner shell transition that lies between 
the intercombination and the forbidden lines of the OVII triplet. We then hypothesized that 



- 21 - 



the A = 13.80 A line be due to NeVIII Ka inner shell absorption. Scahng by wavelength 
ratio from the 'OVII-triplet plus OVI-Ka' system we derived a wavelength for the NeVIII 
Ka of A = 13.78 A, consistent with the observed one within the AA = 20 mA calibration 
uncertainties (see below). Similarly, we tried to estimate the wavelength of the NV Ka inner 
shell transition scaling from the OVII+OVI system, and found X{NVKa) = 29.44 A, in good 
agreement with the HULLAC value of A = 29.46 A. This is significantly different from the 
wavelength of the observed feature at A = 29.54, ruling out the z — NV Ka identification. 
We conclude that at most one of the ten detected hues, at A = 13.80 A, can be tentatively 
identified with a, z — line from NeVIII. 

All other lines must belong to redshifted systems. Nine of the 10 lines can indeed be 
identified as belonging to 2 different intervening absorption systems, at z — (0.011 ± 0.001) 
[cz = (3300 ± 300) km s"^] and z = 0.027 [cz = (8090 ± 300) km 8"^ Table 2]. We focus on 
these two systems in this paper. The strongest line from the absorber at 2; = O.OlliO.OOl, is 
the OVllxa, while for the z = 0.027 ±0.001 absorber the two strongest lines are the NVII^q; 
and the NVlxa- The first five columns of table 2 list, in order of increasing wavelengths, 
measured line wavelengths, equivalent widths (EWs), ion column densities and significances 
for all lines for which a z = identification could not be found . The last three columns 
of Table 2, list the proposed line identifications, rest-frame wavelengths and redshifts. 



3.2.2. Absorption Line Strengths 

We estimated the significance of the absorption lines by two methods: (1) conservatively 
(Table 2, Col. 4), by computing the 3o" error on their EWs, leaving all gaussians parameters 
and the continuum free to vary; (2) by directly integrating the deviations across the line 
profile (Table 2, Col. 5), from the fixed best fitting continuum. Direct integration in all 
cases gives larger significance. The significance of the lines depends critically on the degree 
of accuracy with which the local continuum is known. To assess this dependence for the the 
weakest lines in our spectrum we repeated locally the test on our best fitting continuum 
for all the 'interesting' regions plotted in Fig. 8, after ehminating narrow (50-100 mA) 
intervals of data centered on the z — lines. We let the normalization of our best fitting 



^"Thc' tenth lino is clearly present at an integrated significance of 5a at A = 24.97 A, but could not be 
identified with any known transitions at the rcdshift of any of the three systems. 

^-'^ Included is the line at A = 13.80 A, which we identify here with a NcIX Ka line from the z — 0.011 
system, but that could be at least partly due to NVIII Ka absorption at z = 0. This does not affect our 
conclusion, however, since we consider the EW of this line as a 3a upper limit, due to the super-imposition 
of this feature with a CCD node-boundary. 
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continuum vary locally, and re-fitted the data. In all cases, we obtained slightly lower 
relative normalizations (with value ranging from 0.95% to 0.98%) and unacccptably high 
reduced values, for the given number of d.o.f.. As an example the 21.4-22.5 A interval 
gave a continuum relative normalization of 98.9% and Xri^-'^-f-) ~ 1-4(23), corresponding 
to a probability of exceeding of Pcont{> X^) = 9 %. Including gaussians for the two 
lines at 21.85 A and 22.20 A, and re-fitting the data leaving again the relative continuum 
normalization and the 6 parameters of the gaussians free two vary, gives a relative continuum 
normalization of 1.0 and x^(d.o./.) = 1.0(17), corresponding to Pcont+iines{> X^) — 45 %. 
An F-test between the two models, gives a probability of 10% of exceeding F. All other cases 
gave lower Pcont{> X^) P{> P)- We conclude that our estimate of the local continuum 
level in the regions where the absorption lines are present is statistically acceptable, and so 
confirm the goodness of the estimate of the significances of the lines shown in Table 2 (Col. 
4 and 5). 

The NelX^Q, and OYUkp from the system at z = 0.027, as well as the OVlllKa line 
from the system at 2; = 0.011, although clearly detected in the spectrum, lie close to ACIS 
detector chip node boundaries (Fig. 8a,b). These boundaries may induce an inaccurate 
estimate of the line EWs. We consider the EWs of these lines only as upper-limits. This 
does not affect our main conclusions on the physical state of the two absorbers, as shown in 
sections 7.3 and 7.4.. A 3(7 upper limit for the OVIIIi^-Q from the z — 0.027 system, which is 
not detected in the spectrum, is also listed in Table 2, since it will be used in our diagnostic 
analysis in §6.4. 

Uncertainties in the position of the line centroids are dominated by systematics in the 
LETG wavelength calibrations. Non-linearity in the LETG dispersion relationship produces 
an RMS deviation of predicted vs observed hne positions of 13 mA (Drake, Chung, Kashyap, 
Ratzlaff, 2003), equivalent to 195 km s~^ at 20 A. This uncertainty affects all spectra taken 
with the HRC-LETG configuration, but may be present at a lower extent also in spectra taken 
with the ACIS-LETG configuration. In our combined LETG spectrum, this uncertainty 
easily dominates the few tenths of mA la uncertainty in the gaussian centroid for the 2; ~ 
OVllxa hne (Table 2). We checked the amplitude of the LETG wavelengths accuracy by 
looking for differences between expected and observed relative positions of two strong lines 
from the same ion, in the same system. The largest discrepancy that we find is between the 
position of the OVIIxq and OVIIkp at 2; ~ 0, for which we measure a relative difference of 
AA = 20 mA for a line separation of ~ 3 A. For the absorption hues of the two 2; > systems, 
this is ~ 5 — 50 times larger than the statistical errors on the hues centroid. This calibration 
error fluctuates with wavelength, so that a simple shift is not sufficient to correct for it. We 
have then added this systematic uncertainty to all line centroids wavelength errors in our 
analysis. This translates into average redshift uncertainty ~ Az — 0.001 (i.e. A{cz) ~ 300 



Table 2: Best fitting UV and X-ray absorption line parameters 



Chandra-hFiTG 




X 


EW" 


Column 


Conservative 


Direct Integ. 


Line Id. 


Xrf 


Redshift 


(A) 


(mA) 


lO^'"^ cm-2 


significance 
a 


significance 

a 




A 




13.80 ±0.02 


< 1.5= 


< 1.2= 


[3] 


3.8 


NelX^c. 


13.447 


0.026 ±0.001 


19.11 ±0.02 


< l.S'^ 


< 3.8= 


[3] 


4.0 


OYllKf} 


18.63 


0.026 ±0.001 


19.18 ±0.02 


< 4.r 


< 3.0= 


[3] 


4.9 


OVIIIkc 


18.97'' 


0.011 ±0.001 


19.477 (fixed) 


< 1.8 


< 1.3= 


[3] 


[3] 


OVIIIko 


18.97'' 


0.027 (fixed) 


21.85 ±0.02 


"^■'-'-0.8 


-■■■^-0.3 


3.8 


6.3 


OVII;f„ 


21.602 


0.011 ±0.001 


22.20 ±0.02 


2.2 ±0.8 


0.7 ±0.3 


2.8 


3.8 


OYIlKa 


21.602 


0.028 ±0.001 


25.04 ±0.02 


1.8 ±0.9 


0.8 ±0.4 


2.0 


3.2 


NVIIko 


24.782'' 


0.010 ±0.001 


25.44 ± 0.02 


3.4 ± 1.1 


1.4 ±0.5 


3.1 


5.8 


NVIIk„ 


24.782'' 


0.027 ±0.001 


29.54 ±0.02 


3.6 ± 1.2 


0.7 ±0.2 


3.0 


4.6 


NVIko 


28.787 


0.026 ±0.001 


34.69 ±0.02 


2.4 ± 1.3 


0.5 ±0.3 


1.8 


3.3 


CVI^„ 


33.736'' 


0.028 ±0.001 


FUSE 




1042.41 (fixed) 


< 18= 


< 0.014= 


[3] 




OVl2«^2p 


1031.926 


0.01 (fixed) 


1043.28 (fixed) 


< 21 


< 0.016 


[3] 




OVl2s^2p 


1031.926 


0.011 (fixed) 


1059.79 (fixed) 


< 19^ 


< 0.014= 


[3] 




OVl2.^2p 


1031.926 


0.027 (fixed) 


HST-GHRS 




1228.02 ±0.01 


92.0 ±3.3 


0.0180 ±0.0006 


27.9 


97.2 


Hlij/a 


1215.67'' 


0.010160 ±0.000007 


1229.04 (fixed) 


< 26= 


< 0.0047= 


[3] 






1215.67'' 


0.011 (fixed) 


1248.85 (fixed) 


< 47= 


< 0.0085= 


[3] 






1215.67'' 


0.027 (fixed) 



Absorption fine EW are formaUy negative; here we drop the sign. '' Central wavelength of the two fines of the doubled, computed by 
weighting their positions by their oscillator strengths. = 3a upper limit. 
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km s-^) at 20 A. 

Figures 10 and 11 show all the detected metal hnes from our two intervening systems in 
velocity space as residuals (in a) between the best fit continuum and the data. The vertical 
dashed hne across the panels is the average redshift of the system. Because of the systematic 
LETG wavelength errors the hne wavelengths have been shifted to the mean z of the system. 
These difference are always within the estimated systematics in the wavelengths calibration 
(i.e. AA < 20 mA). 

3.2.3. Statistical Significance of the Two Absorption Systems 

We estimated the significance of the two systems z — 0.011 and z — 0.027 by 
comparing the best-fit of simpler models that do not include the metal absorption lines 
from these systems with that of more complex models that include them. The F-tests 
were run over the following two portions of the spectrum, containing a similar number 
(100 and 130) of resolution elements: AX{z = 0.011) = [21,26] A and AX{z = 0.027) = 
[21.5,23] U[24.5,26] U[27.5,30] U[34,35] A. OVllRa and NVllxa lines were included in the 
complex model for the z = 0.011 system, and OVIIj^q, NVllxa, NVIi^Q, and CVIj^q were 
used for the z — 0.027 system. The fits with the complex models were performed freezing 
the width of the lines to the instrument resolution and fixing their relative positions and 
normalizations to their expected relative position based on the proposed identifications, 
and their measured relative normalization. Only the redshift of the system and the overall 
normalization were left free to vary in the fits. The results of these tests are presented in 
Table 3, in which we list both the probability of exceeding F (Col. 3) and the associated 
significance in a (Col. 4). The last column of Table 3, lists the sum of the lines significances 
for the two systems (Col. 4 of Table 2), computed only over those lines detected individually 
with a significance > 2a. 

Table 3: Significance of the Absorption Systems 



Redshift 


> 2(7 Ion 


Probability of 


Significance 


Sum of > 2a lines 




Detections 


Exceeding F 


in a 


significances (in a) 


0.011 ±0.001 


OVII. NVII 


5.9 X 10-4 


3.5 


5.8 


0.027 ±0.001 


OVIL XVIL XM 


1.3 X U)-'' 


4.9 


8.9 
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Fig. 10. — Residuals in a, in velocity space, for three absorption lines of the cz = 3300 ± 300 km s"^ 
system. 
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Fig. 11. — Residuals in cr, in velocity space, for six absorption lines of the cz = 8090 ±300 km s ^ system. 
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4. The HST and FUSE Spectra of Mkn 421 
4.1. HST-GHRS 

Mkn 421 was also observed with the Hubble Space Telescope (HST) GHRS on 1995 
February 1. An intervening HI Lya system at cz — (3046 ± 2) km s~^ was discovered and 
reported by ShuU et al. (1996). We retrieved this GHRS spectrum of Mkn 421 from the 
public HST archive, and re-analyzed the data. Wc fitted the 1220-1250 A HST spectrum 
with a combination of broad continuum components and three narrow absorption gaussians 
to model the intervening HI Lja and two strong Galactic Si lines visible in the spectrum 
(Figure 12, Table 2). Our best fitting parameters for the HI Lya line (we assume the observed 
b = iOtl km s-i - errors are at 3a - to estimate the HI column) are consistent, within the 
uncertainties, with those reported by Penton et al. (2000a). The redshift of the HI Lyct fine 
(Tab. 2) is consistent with that of the z — 0.011 X-ray absorption system within the large 
systematic X-ray uncertainties {Az — 0.001). We also computed the 3a HI column upper 
limit of a putative HI Lya at the average redshifts of the two z > X-ray systems (Tab. 2), 
and with b{z = 0.011) = 130 km s"^ and b{z = 0.027) = 185 km s'^ (i.e. T{z = 0.011) = 10^ 
and T{z = 0.027) = 2 x 10^ K: see §7.4). We found Nhi{z = 0.011) < 4.7 x 10^^ ^^-2 
Nhi{z = 0.027) < 8.5 x lO^^ ^^-2 (j^^^^ 2). 



4.2. FUSE 

Following our TOO request Mkn 421 was observed by the Far Ultraviolet Spectrometer 
Explorer (FUSE) on 2003, January 19-21 (MJD = 52658-52660), with a total exposure time 
of 62 ks, based on a trigger from another period of strong X-ray activity as recorded by the 
RXTE-ASM (Figure 1). 

We searched the FUSE spectrum of Mkn 421 for OVl2s^2p absorption at or close to the 
redshifts of our two z > X-ray absorbers (Figure 13a,b). Wc fitted the two portions of the 
FUSE spectra between 1040 and 1051 A (Fig. 13a) and 1056 and 1070 A (Fig. 13b) with 
a combination of a polynomial continuum and a number of absorption gaussians to model 
the absorption features (mainly Galactic H2 lines) visible in these two spectral regions. In 
particular, we fitted 2 H2 lines and one Arl fine between 1040 and 1051 A (Figure 13a), and 
one H2 and one Arl line between 1056 and 1070 A (Figure 13b). We then added to our 
models two absorption gaussians with same FWHM, with their relative positions forced to 
the separation of the OVl2s_>2p doublet lines (A(0Vli)/A(0Vl2) = 0.9945), and with their 
relative amplitudes forced to be 2, the relative amplitude of unsaturated OVl2s-+2p doublet 
lines (i.e. Ampl(0VIi)/Ampl(0Vl2) — 2). For the X-ray system at z=0.011, the strongest 
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Fig. 12. — HST-GHRS spectrum of Mkn 421 between 1223 and 1256 A. The best fitting continuum plus 
absorption line model is shown as a solid line. The two vertical arrows mark the positions of the HI Lya at 
the redshifts of the two z > X-ray absorbers. The horizontal error bars on these arrow show the ±300 km 
s~^ systematic uncertainty of the LETG. 
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Fig. 13. — Two portions of the FUSE spectrum of Mkn 421 in the 1040-1051 A (upper panel) and 1056-1070 
A (lower panel) ranges, along with their best fitting continuum plus absorption line models. The shaded area 
in Fig. 13a, indicates a detector flaw, due to a dead spot in the detector (Savage et al., 2005). The horizontal 
error bars in the two panels show the position of the expected lines of the OVI doublet at z = 0.011 (Fig. 
13a) and z = 0.027 (Fig. 13b) ±300 km s'^ systematic uncertainty of the LETG. 
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line of the OVl2s^2p doublet would fall at A = 1043.28 A, where a detector artifact due to 
a dead spot (Savage et al., 2005; Figure 13a) hampers its detection. The position of the 
second line of the putative doublet would fall at A = (1049.03 ± 0.05) (Figure 13a) where 
a moderate deficit in the continuum, compared to the best fit model, is visible. However, 
this region of the spectrum is potentially contaminated by Galactic H2 absorption (4-0 R(0) 
at A = 1049.367), so making the identification of a putative A = 1037.617 A OVI line at 
z — (0.011 ± 0.001) unreliable. Unfortunately, we cannot use the FUSE spectrum of Mkn 
421 to measure OVI absorption z — (0.011 ± 0.001) and therefore we only measure a 3(T 
upper limit on the OVI column of ^ovi{z = 0.011) < 1.6 x 10^^ cm-^ (Table 2). No OVI 
lines are visible at the redshift of the HI Lya absorber, which gives a 3(7 upper limit on 
the column density of OVI at that redshift of Nov/(-2 = 0.01) < 1.4 x lO^''^ cm'^ (Tabic 2). 
Similarly, we do not detect any OVI counterpart for the X-ray system a,t z — 0.027, down to 
a 3a sensitivity of Nqvi < 1-4 x 10^^ ^,^^-2 (^pigure 13b, Table 2). 



5. Summciry of Absorption System Properties 

5.1. The z — 0.011 Absorption System 

The z = 0.011 X-ray system is mainly detected through OVII and NVll absorption. No 
NelX, NVI or CVI absorption lines are detected at this redshift, while only a weak 3a upper 
limit can be set on OVIII absorption, due to the presence of a strong instrumental feature. 
Based on the detected lines we measure an average redshift of the system of 2; = O.OlliO.OOl, 
or cz = (3300 ± 300) km s'^ (Table 1). Due to an instrumental flaw affecting the FUSE 
spectrum of Mkn 421 between A = 1043 A and A = 1043.8 A (Figure 13a), only a 3a upper 
limit can be set on the column of OVI ai z = 0.011. Finally, a strong HI Lya is detected at 
a consistent redshift in the HST-GHRS spectrum of Mkn 421 (ShuU et al., 1996). 



5.2. The z = 0.027 Absorption System 

The z — 0.027 system has four identified detections of resonant absorption transitions 
from as many different ions (CVI, NVI, NVII and OVII), plus two uncertain EWs (because 
of the presence of two instrumental features) detections of NelX and the OVUkis, which we 
treat as 3a upper limits. Based on the four detected lines we measure an average redshift of 
the system of 2; = 0.027 ±0.001. No HI Lya or OVl2s^2p absorption is detected in the GHRS 
and FUSE spectra of Mkn 421 at a redshifts oi z = 0.027. Assuming a temperature of ~ 10^ 
K and pure thermal broadening of these lines (plus stretching due to the Hubble flow for an 
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overdensity extending 2 Mpc along our line of sight) , we estimated 3a upper limits of 47 mA 
and 18 mA respectively for these transitions, corresponding to Nhi < 8.5 x 10^^ cm~^ and 
Novi < 1.4 x lO^^ cm'^ (Table 2). 

6. Local Galaxy Environment of the X-ray Absorbers 

To evaluate the density of the galaxy environment in which the X-ray absorbers lie, 
we searched a number of galaxy catalogs sensitive down to = 15.5: the CfA Redshift 

Catalog (v. June 1995, CfA95 ), the CfA-North-Galactic-Pole+36 Catalog (CfA-NGP+36, 
Huchra et al., 1995, ApJS, 99, 391) and the Updated Zwicky Catalog (UZC99, Falco ct al., 
1999, PASP, 111, 438). The CfA-NGP+36 completes, in redshift identifications, the CfA95 
catalog in the declination range 32.5°-38.5° and right ascension range 8^ < « < 17^* (which 
includes the position of Mkn 421), and down to nib = 15.5. The UZC99 catalog corrects 
some of the inadequacies of the previous CfA catalogs, which were all based on the Zwicky 
Catalog of galaxies and relied on heterogeneous set of galaxy positions and redshifts. 

For the z = 0.011 system, we found 4 galaxies (2 of which forming a close system, in 
interaction) in a comoving cylinder with a base-radius of 5 Mpc, and half-depth of 500 km/s 
(in the z dimension: 2800 < cz < 3800 km s-^ Figure 14; see also Penton et al, 2000b) 
. All these galaxies have A{cz) < 140 km s~^ from the z — 0.011 system (within the absorber 
redshift uncertainty), and he between 4.4 and 5.1 Mpc from it. At the absorber distance, the 
limit magnitude of the catalogs, translates into an absolute magnitude limit of Mf, = —17.9, 
i.e. 2.75 mag fainter than M*j^{SDSS) (Poli et al., 2003 ). This implies a density of 
galaxies at ;z = 0.011, down to Mj, = -17.9, of 0^^,°-°"(Mb > -17.9) = 3.6 x 10"^ We used 
the transformed g* to B band SDSS local Luminosity Function of galaxies (see Fig. 1 in 
Poli et al., 2003), to compute the local density of galaxies down to this magnitude limit, and 
found (psDssi^b < —17.9) = 0.017 Mpc~^. This gives a relative density of galaxies in our 



^2http://cdsweb.u-strasbg.fr/viz-bin/Cat?VII/193 

^•^This depth correspond to a physical size of the cyhnder in the z direction, of ~ 14 Mpc, if pecuhar 
galaxy motions are neglected. 

^^The linear sizes of this volume, correspond to the typical sizes expected for the predicted intergalactic 
filaments connecting virialized structures in the local Universe (i.e. Cen & Ostriker, 1990, Hellsten et al., 

1998). 

^^Poli ct al. (2003) derive M*g{SDSS) ~ —20.65 from the g magnitude local Luminosity Function of 
galaxies from the Sloan Digital Sky Survey (SDSS, Blanton et al., 2001) adopting the following transformation 
Mb = g * -1-0.32 and scaling magnitude and normalization to Hq = 70 km s~^ Mpc~^. 
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region of Universe centered on the z = 0.011 system, of 0^^^° °^^(M6 > -\1 .^)l(^sDSs{Mb < 
—17.9) = 0.2. The g* band SDSS luminosity function is an average over 4684 galaxies in a 
volume ~ 5.1 x lO^/ig f Mpc^ (Blanton et al., 2001) and so is a good mean sample of the 
local Universe. We conclude that this system lies in an under dense volume of Universe. 

The z = 0.027 system lies ~ 13 Mpc from Mkn 421. The CfA95, CfA-NGP+36 and 
UZC99 catalogs contain a total of 8 different galaxies down to = 15.5 in a comoving 
cyhnder with a base-radius of 5 Mpc and half-depth of 500 km/s (in the z dimension: 
7600 < cz < 8600; Figure 15)^^ is r^^ggg galaxies lie between 3.3 and 8.2 Mpc from the 
absorber. At the absorber distance, the limiting magnitude of the catalogs, translates into 
an absolute magnitude limit of = —19.8, or 0.85 mag fainter than M*^{SDSS) (Poli et 
al., 2003). The observed density of galaxies at z = 0.027, down to this magnitude limit, is 
then 0^^° °^''(Mft < —19.8) = 7.3 x 10"^ Mpc~^. The local density of galaxies for the same 
magnitude hmits, from the B band SDSS local luminosity function of galaxies (Poli et al., 
2003), is also ^sossiM^ < -19.8) = 7.3 x 10~^ Mpc"^, which implies a relative density of 
galaxies of 0^;;° °2^(M6 < -19.8) /(psnssiMb < -19.8) = 1. The z = 0.027 system seems to 
lie in an average density region. 



7. Discussion 

7.1. Line Saturation and EW Ratio Diagnostics 

In neither X-ray absorption system do we detect more than one line from the same ion. 
So we can not use the standard EW ratios technique to estimate line saturation directly. Wc 
instead evaluated line-saturation by using the standard alternative Curvc-of-Growth (CoG) 
technique. We allowed the Doppler parameter to vary between 6 = 50 and b = 400 km s^^. 
Since all the lines are unresolved we can place a 2a limit of 6 < 480 km s~^. The largest 3a 
determinations of our hue EWs, for both systems, fell well within the hnear branch of the 
corresponding CoGs even for 6 = 50 km s~^, and far from saturation. 

On the linear branch the ratio of the EWs of two lines from two ions X* and of 
the two elements X and Y, is directly proportional to the ratio between the corresponding 
relative ion abundances ix^/^YJ- 

_ Wx^jmA) n{Y^) y fy. Ay 

. In eq. 1 fyj / fx^ is the oscillator strength ratio for the given transition, and Ay/ Ax is 
the relative abundance of element Y compared to X. Two or more of these ratios from ions 
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Fig. 14. — Aitoff plot of a region of sky, in Galactic coordinate around the z = 0.011 X-ray absorber. 
Four galaxies are found in a cylinder centered on the absorber, with a base radius of 5 Mpc and half-depth 
of A{cz) = 500 km s~^ (~ 7 Mpc), down to a magnitude limit of mt = 15.5. Squares and circles indicate 
blueshifted and redshifted galaxies respectively, compared to the absorber redshift. The star marks the 
position of the absorber. 
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Fig. 15. — Aitoff plot of a region of sky, in Galactic coordinate around the X-ray absorber at z — 0.027. 
Eight galaxies are visible in a cylinder centered on the absorber, with a base radius of 5 Mpc and half-depth of 
A(cz) = 500 km s~^ 7 Mpc), down to TOf, = 15.5. Squares and circles indicate blueshifted and redshifted 
galaxies respectively, compared to the absorber redshift. The star marks the position of the absorber. 
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of the same element, then, uniquely define the ionization balance of the absorber, and so 
allow us to derive the physical conditions in the gas, as well as the ionization mechanism at 
work. The ratio of the EWs of lines from ions of different elements also depends on the ion 
abundance and so constrains the chemical composition of the gas. 

We used Cloudy (vs 90.04, Ferland et al., 1997) to build grids of hybrid, coUisional 

ionization plus photoionization models, as a function of the electron temperature, to compare 
predicted ion relative abundance ratios with observed EW ratios. These models reduce to 
pure coUisional models for electron density of the gas >> 10~® — 10~^ cm~^. 



7.2. Modeling the Photoionization Contribution 

Following Nicastro et al. (2002), we used the diffuse UV-to- X-ray background (XRB, 

hereinafter) as the photoionizing source. Wc also included the possible 'proximity effect' (c.f. 
HI Lya, Bechtold, 1994) of photoionization by the blazar Mkn 421. A 100 mCrab flux for 
Mkn 421 implies beamed X-ray luminosities along our line of sight, of the order of LMkn 421 ~ 
]^q45 gj,g g-i^ Such luminosity can photoionize low density gas along our line of sight to 
Mpc-scale distances. Assuming the observed UV-to-X-ray Spectral Energy Distribution 
(SED) of Mkn 421, the ionization parameters , at the illuminated face of the clouds, are 
Umuji 421 — 1.7 X 10~^(5~^ at 2; = 0.011 and Umuu 421 — 0.8(5"-^ at 2; = 0.027 (we have used 
Jib = 0.04 and h — 0.7). The recombination time for OVII in photoionized gas with out-of- 
equilibrium temperature of T = 10^ K and overdensity 5^500 is tree > 2.4 x lO^^^g years 
(e.g. Nicastro et al., 1999, and references therein), and increases/decreases for lighter /heavier 
elements. The inverse, equilibrium 'photoionization' time is typically more than one order 
of magnitude shorter, and depends on the variability timescale of the ionizing source. Given 
the frequency of the 100-mCrab type bursts recorded from Mkn 421 (~ 1.5/year in the past 
7 years), and their typical duration of ~ 3 weeks, we must conclude that any intervening 
IGM condensation within ~ 10 — 20 Mpc of this blazar with overdensities 5:^500 must be 
kept in a relatively high ionization state due to this proximity effect. Moreover, given the 
different times with which different ions react to ionizing flux changes the overall ionization 
balance in such a gas could be significantly different from the one expected at equilibrium. 
The photoionization contribution due to the diffuse XRB can be parameterized, for both 
absorbers, by the ionization parameter Uxrb ^ 0.15 x (again, assuming (0^/1^/0.02) = 1, 



"^^U is here defined as the ratio between the radiation field ionizing photon density at the distance of 
the gas cloud, and the electron density of the gas, and is computed assuming constant density through the 
cloud, and a plane parallel geometry (i.e. thickness of the cloud along the line of sight, much smaller than 
the distance of the cloud from the ionizing source). 
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and neglecting the small redshift difference between the two systems). By comparing the 
proximity effect with the ionization induced by the UV-X-ray background for the two systems 
we conclude that photoionization by Mkn 421 is completely negligible for the filament at 
z = 0.011 , while it dominates the overall photoionization contribution by a factor of ~ 5 
for the filament at z = 0.027. 



7.3. Self-Consistent Nh versus T Solutions 

For both the z = 0.011 and z = 0.027 absorption systems we compared the measure- 
ments (X-rays) or 3a upper limits (UV, FUV and X-rays) on the EW ratios (§7.1), with 
the predicted ion relative abundance ratios from our grid of hybrid collisional-ionization plus 
photoionization models (§7.2), to constrain maximal intervals of temperatures defining com- 
mon solutions (e.g. Nicastro et al., 2002). The existence of a common range of temperatures 
for all ion ratios, however, still does not guarantee that a range of self-consistent equivalent 
hydrogen column density versus temperature solutions exist. To search for these solutions we 
used the each measured (or 3a upper limit) equivalent width Wx for a given ionic transitions 
to estimate a temperature (and metallicity) dependent equivalent hydrogen column density 
^h{T, [X/H]) (where [X/H] is the logarithm of the metallicity of the generic element X, 
compared to solar). The temperature dependence is introduced by the dependence of Nh 
upon the relative ion abundance ^i, in the versus Wx relationship (see e.g. eq. (1)). 
This way we were able to define two ranges of self-consistent Nh versus T solutions for the 
two systems, and constrain their metallicity (§7.4 and 7.5). 

7.4. Physical State of the Absorber at 2; = 0.011 

The only X-ray absorption line of the z = 0.011 system individually detected at a 
significance larger than 3a is the OVllKa fine at 21.85 A (Table 2). So we use this ion as 
the reference ion for the ionization balance diagnostics. 

The predicted ratios for several ions compared to OVII, for hybrid collisional-ionization 
plus photoionization gas, are shown as thin solid curves in Figure 16, plotted against the 
gas temperature: (a) OVI/OVIl (black curve), (b) OVIlI/OVll (red curve), (c) NVII/OVII 
(blue curve), and (d) HI/OVII (green curve). The predicted ratios are for a baryon density 
of Ub — 10~^ cm~^ {S — 50). The thick segments superimposed on these curves are our ±3a 
allowed ranges for the corresponding ratios. We assume [0/H] = -1 and [N/0] = 0. Lower 
values of [0/H] would further reduce the allowed range of temperatures set by the HI/OVII 
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interval (green segment on the green curve), while larger [0/H] would broaden the interval. 

First we note that the observed width of the HI Lya at A = 1228.02 ± 0.01 (Figure 12, 
Table 2) sets a 3a upper limit on the temperature of the HI absorber (dashed vertical hne 
in Figure 16) of Thi — 1.2 x 10^ K. This is incompatible with the, metaUicity independent, 
lower hmit of Tpuv-x-ray — 5.3 x 10^ K on the temperature of the FUV- X-ray absorber, 
imposed by the measured 3cr lower boundary on the OVI/OVII ratio (Fig. 16, left end of the 
OVI/OVII black segment). Decreasing the overdensity of the FUV- X-ray absorber from the 
adopted value of 5 = 50 to, say, S = 10, would slightly flatten the OVl/OVll curve at low 
temperatures, extending the OVl/OVll lower end interval to Tpuv-x-ray ~ 3.9 x 10^ K. This 
reduces the inconsistency between the temperatures of the HI and FUV-X-ray absorbers, but 
they remain mutually incompatible. We conclude that if the HI and FUV-X-ray systems 
are related, they must co-exist in a multiphase IGM. In the following we use the HI Lya 
3(7 upper limit at A = 1229.04 k {z — 0.011, Table 2) as a constraint for condition in the 
FUV-X-ray absorbing gas. 

Analogously to the OVI/OVII ratio (which sets a minimal temperature of T'^uv -x-ray = 
5.3 X 10^ K), the measured OVIII/OVII 3a upper limit sets a metaUicity independent max- 
imal temperature of the absorber at Tp^y_x-ray ^ 3.3 x 10^ K (Fig. 16, right end of the 
OVIII/OVII red segment). 

On the contrary the measured {3a) HI/OVII and NVII/OVII ratios depend upon metal- 
hcity. In Figure 16 the thick green and blue segments refer to Solar and O.lx Solar metaUicity 
respectively. These measurements define broad ranges of temperatures that both intersect 
the metaUicity independent maximal range defined by the OVI/OVII and OVIII/OVII ratios. 

Additionally, the 3a upper hmit on the NVII/OVIII ratio constrains the [0/N] relative 
abundance to be [0/N]> —2.05, consistent with solar . Similarly, and more importantly, 
we can derive a lower limit on the oxygen metaUicity ratio. The non-detection of HI Lya 
together with the detection of OVll Ka and the minimum amount of relative HI to OVII 
abundance predicted by the model, set a stringent lower limit on the oxygen to hydrogen 
metaUicity ratio, of [0/H]> -1.46. 

The two arrows in Figure 16 shows the extremes of the maximal range of temperatures 
of the z = 0.011 FUV-X-ray absorber: Tpuv-x-ray G [0.53 - 3.3] x 10^ K. Coupling this 
maximal temperature range with the estimates on Nh obtained from each available ion 



'From eq. 1, the relative metaUicity between two elements Y and X depends on the inverse relative ion 
abundance ratio ix^/^Yi^ and on the direct EW ratio Wyj/Wxi- For a Wxi 3(j upper limit, this relative 
metaUicity has to be always larger than the 3a lower limit on the ratio Wyj /Wx* , multiplied by the minimum 
allowed value of the curve ^x' /^n ■ 
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Fig. 16. — Predicted ratios for several ions compared to OVII, for hybrid coUisional-ionization plus pho- 
toionization gas (thin sohd curves), plotted against the gas temperature: (a) OVI/OVII (black curve), (b) 
OVIII/OVII (red curve), (c) NVII/OVII (blue curve), and (d) HI/OVII (green curve). The predicted ratios 
are for a baryon density of rib = 10~^ cm~^ {S = 50). The thick segments superimposed on these curves are 
our ±3(T allowed ranges for the corresponding ratios for the z — 0.011 system. We assume [0/H] = -1 and 
[N/0] = 0. The dashed vertical line sets a 3a upper limit on the temperature of the HI absorber, from the 
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(see §7.3), we find that a common range of Nh versus Tpuv -x-ray solutions exists in the 
temperature interval Tg G [0.6-2.5], and is given by: Nh(^ = 0.011) = [(1.7±0.7)T6- (0.1± 
0.6)] X 10^^ X 10"['^/^l-i cm~^ (where Tg is the temperature of the absorber in units of 10^ 
K). This implies a thickness of the absorber along our line of sight D{z = 0.011) =^H/nb = 
[(0.55 ± 0.23)T6 - (0.03 ± 0.19)](n_5b)-il0-[^^/^l-i Mpc, assuming constant baryon density. 
In this equation (n_5;,) is the baryon density in units of 10^^ cm~^. 

We notice here that due to the weakness of the OVlll/OVll and NVll/OVll ratio 
constraints, our ionization balance solutions are rather insensitive to even factor 10 differences 
in the baryon density. Higher signal to noise X-ray data are needed to constrain both the 
OVIII and the NVII columns to better than 30-40 % at 2a, and so set limits on the value of 
the baryon density. This is obtainable with 2.5-3 times the number of counts we currently 
have at the wavelengths of the OVIII and NVII Ka transitions. 



7.5. Physical State of the Absorber at z ^ 0.027 

Two X-rays lines of the z = 0.027 system are detected at a conservative significance 
> 3a: the NVIIj^^ at A = 25.44 A and the NVI^a at A = 29.54 A (Table 1). In our 
ionization balance diagnostics analysis we use NVII as the reference ion as it has slightly 
higher significance. 

As for the z = 0.011 system we explore the hybrid coUisional ionization plus photoion- 
ization case for the z — 0.027 system (Figure 17). Seven EWs ratio constraints from 8 
different ions are now available from the UV-FUV- X-ray determinations or 3a upper limits 
(Figure 17). 

First we note that the NVI/NVll constraint (blue interval in Figure 17) sets a tight and, 
most importantly, metallicity independent lower limit on the temperature of the absorber, 
of logT> 5.98 (vertical black solid fine and blue arrow in Fig. 17). Based on the simple 
equilibrium models we adopt, however, a maximal temperature range of solutions can be 
found only for rather large values of the N/0 ratio, compared to Solar. Solid intervals on 
the 0/NVII curves in Fig. 17 are for [N/0] = 0, while broader dashed intervals are for 
[N/0] = 0.65. This is the minimum N/0 overabundance, compared to the solar value, 
that reconciles the large observed absorption by He-like and H-like N with the lack, or 
weakness, of OVl, OVll and OVlll absorption. In particular, from the high temperature 
side (logT>6.4 K), [N/0]> 0.65 allows the measured OVI/NVII and OVII/NVII upper limit 
ratios to exceed the predicted OVII/NVII local maxima at the knee temperature of T ~ 10^-^ 
K, where the curve changes its slope and becomes much flatter. On the other hand, from 
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the low-temperatures side (logT< 6 K), only for [N/0] ^0.65 can the OVIII/NVII constraint 
extend up to the logT> 6 K region, and so become compatible with both the OVI/NVII 
and the OVII/NVII constraints. A large N abundance is also consistent with the somewhat 
narrow range of allowed temperature that is derived from the Hl/NVll measurement. The 
red interval on the theoretical Hl/NVII curve is for [N/H] = -0.35, which together with [N/0] 
— 0.65, implies [0/H] = -1. At this relatively large (but still sub-Solar) [N/H] metallicity 
ratio, the temperature of the system is constrained between 5.97 < logT < 6.49 K. Larger 
or lower [N/H] values would broaden or narrow this interval. The minimum N/H metallicity 
ratio for this system is [N/H]> -0.67 (i.e. [0/H]> -1.32, for [N/0] = 0.65) at logT = 
6.2 K. Finally the N/C ratio is consistent with the Solar value while only a moderate N/Ne 
overabundance of [N/Ne] = 0.08 is required (so maybe indicating that also C and Ne are 
particularly abundant in this absorbed, compared to H). 

We stress that the above limits on the metaUicity ratios depend critically on modeling, 
and so on (a) the density of the absorber, and (b) whether or not equilibrium photoionization 
(acting on gas already collisionally ionized through shocks) applies, and applies to all ele- 
ments (we recall that recombination times are different for different elements). The predicted 
relative abundance ion ratios plotted in Fig. 17 are for a baryon density of rih = 10^^ cm~^ 
{S = 47). Larger densities would generally increase all predicted ratios, so making the N 
overabundance problem even more serious. Lower densities, however, would generally lower 
the OVI/NVII and OVII/NVII knees at logT~ 6.4 K, but would increase the OVIII/NVII 
ratio at low temperature: a,t rib — 2 x 10~^ cm~^ (S ~ 10), a [N/0] = solution can be 
found for all ions, except OVIII, with a temperature of logT~ 6.2 K. Finally, the gas could 
be far from equilibrium, in which case no conclusion could be drawn of course on metallicity 
and temperature of the absorber based on our comparison with hybrid equilibrium models. 

For [N/H] — —0.35 and [N/O] > 0.65 the maximal temperature range for which Nh 
versus T are in principle allowed is set, on the lower side by the metallicity independent 
NVI/NVII ratio, and on the upper side by the metallicity dependent HI/NVII ratio: Tq E 
[0.95 — 3.1]. Coupling this maximal temperature range with the estimates on N// obtained 
from each available ion (see §7.3), we find that, for [N/O] = 0.65, a common range of Nfj 
versus T solutions exists in the narrow (but, again, metallicity dependent) temperature 
interval Tg G [1.2 - 1.6], and is given by: Nniz = 0.027) = [(1.67± 0.76)T6 - (0.71 T 1-15)] x 
10^^ X lO~i\-0/H]-i+iN/o]o.66} cm"^. This implies a thickness of the absorber along our line of 
sight D(z = 0.027) ^Nn/nb = [(0.53 ± 0.24)^6 + (0.23 ^ 0.37)](n_56)-H0-{[^/-^l-i+[^/^l° «5} 
Mpc, assuming constant baryon density. 
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Fig. 17. — Same as Fig. 16, for the z — 0.027 system. Predicted ratios for 7 different ions are compared 
NVII, for a baryon density of = 10"^ cm^^ ^ 50). We assume [0/H] = -1 and [N/0] = (solid 
interval) and 0.65 (dashed intervals). The NVI/NVII constraint (blue interval) sets a tight and metallicity 
independent lower limit on the temperature of the absorber, of logT> 5.98 (vertical black solid line and blue 
arrow) . 
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7.6. Two WHIM Filaments at z ^ 0.011 and z = 0.027 

The z = 0.011 and z = 0.027 absorbers that we detect along the hne of sight to Mkn 421 
have temperature ranges in good agreement with the theoretical prediction for the WHIM. 
The WHIM temperature distribution peaks at ~ 4 x 10^ K and is highly skewed toward 
lower temperature, with T{a_)/Tpeak ~ 27 and T{a^)/Tpf.ak ~ 5 (Dave et al., 2001: see their 
Fig. 5). The hmits that we set on the equivalent H columns of these two absorbers (^niz — 
0.011) = [(1.7±0.7)r6-(0.1±0.6)] X 10^9 X iq-IO/hu ^.^-2 ^^^^ ^ o.027) = [(1.67 ± 
0.76)T6 - (0.71 T 1-15)] x 10^*^ x IQ-iio/Hl-i+lN/oWa^} cm'^) are also in excellent agreement 
with WHIM predictions (e.g. Hellsten et al., 1998, Dave et al, 2001). Predictions for WHIM 
metallicities are still uncertain, mostly due to the lack of, or partial, non self-consistent, 
treatment of galaxy-IGM feedback mechanisms in current hydrodynamical simulations for 
the formation of large scale structure. However, values of about [0/H] = -1 or even larger 
are plausible at very low redshift, consistent with the [0/H] ~ —1 we have found. Moreover, 
the number of detected OVII systems is consistent, within the la errors, with the predicted 
number of OVII WHIM filaments along a random hne of sight, down to the lowest detected 
OVII column of ^ovii — 7x 10^^ cm~^, and up to the distance of Mkn 421 (see next section). 

We therefore identify the two absorption systems at 2; = 0.011 and z — 0.027 with two 
intervening WHIM filaments, with roughly same column and extent along our line sight. 
These are the first > 3.5(7 (conservative, or > 5.8(7 - sum of lines significance) detections of 
WHIM filaments in the X-rays. 

7.6.1. Ruling Out Blazar Outflows 

We should also consider the possibility that the two absorbers are due to highly ionized 
material outflowing from either the nucleus of Mkn 421 or the blazar 's host galaxy at veloc- 
ities of ~ 5600 km s~^ and ~ 900 km s~^ respectively. Even the lower of these two velocities 
is extreme for ISM clouds, when compared with the typical range of velocities observed in 
High Velocity Clouds (either cold, e.g. Blitz et al., 1999, or warm-hot, e.g. Sembach et 
al., 2003) in our Galaxy. Few hundreds to 1000-2000 km s~^ photoionized outflows are, 
instead, common in low redshift Seyfert Is (e.g. Kriss, 2002), but have not been seen in 
blazars. Early reports of X-ray absorbers in blazars (Kruper & Canizares, 1982; Madejski 
et al. 1991) have not been confirmed by later, higher spectral resolution, observations (e.g. 
Nicastro et al., 2002, Fang et al., 2003, Cagnoni et al., 2004). Moreover all photoionized 
outflows in Seyfert Is are detected in absorption both in the UV (through CIV and OVI) 
and X-rays (through hundreds of transitions from ionized metals, e.g. Krongold et al., 2003, 
and references therein) , while in neither of our two cases we do detect associated OVI. 
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In the following we investigate this possibility further. If the absorber was part of the 
blazar environment, assuming a range of density of 10^ — 10^^ cm~^ (typical of different 
gaseous components in an AGN environment, from the parsec-scale dusty torus, to the sub- 
parsec scale Broad Emission Line Region clouds - BELRs - or Warm Absorbers - WAs) and 
distances from 1 pc (the size of the torus for an object with the non-beamed luminosity 
of Mkn 421, i.e. L ~ 10*^^ erg s~^) down to 0.5-1 hght day (the size of the BELR and, 
possibly, the WA) the ionization parameter U at the illuminated face of the cloud, during 
100 mCrab outbursts of the source, would reach value ranging between Utorus — 30 and 
Ubelr-wa — 0.3 — 30, putting most of the oxygen in the form of OVII and OVIII (for the 
BELR cloud) or OVIII and OIX (for the WA or the torus; see e.g. Figure 1 in Nicastro et 
al., 1999). In the spectrum of Mkn 421 we only detect OVII, and measure upper limits for 
OVIII implying f/^ 1, so absorbing material from either the molecular torus or a WA in the 
nuclear environment of Mkn 421 along our line of sight is unlikely to produce the observed 
features. A particularly dense BELR cloud, instead, could. However, typical equivalent H 
column density for the BELR clouds are in the range 10^^ — 10^^ cm^'^ (e.g. Blandford, 
et al., 1990), and metaUicities are typically at least solar (e.g. Korista et al., 1996). This 
would imply OVII columns of > 2 x 10^^ cm~^ (assuming an OVII fraction of fovii — 0-2 at 
Ubelr — 0.3, see Figure 1 of Nicastro et al., 1999), more than 3 orders of magnitude larger 
than the value we actually measure. 

The same argument can be used to estimate the ionization degree of an extremely high 
velocity ISM cloud, located at kpc-scale from the blazar's host galaxy center (and along the 
line of sight to Mkn 421). Assuming typical ISM densities of ;^ 1 cm~^ we obtain C/^30 during 
outburst phases, which implies virtually no OVII, with foix — 0.95 and foviii — 0.05 (e.g. 
Nicastro et al., 1999). 

Finally we also looked for variabihty of the strongest hnes of these two absorption 
systems in the two separate ACIS-LETG and HRC-LETG spectra of Mkn 421, taken about 
6 months apart and in different flux states. Variability of the absorption lines found here 
would denote an AGN origin. The OVII absorber at z = 0.011 is clearly present, even if 
at lower significance, in both the ACIS-LETG and the HRC-LETG spectra and with same 
W ^ 3 mA, If the absorber was intrinsic to the nucleus, and if its density was ng^W^ 
cm~^ (Nicastro et al., 1999) the change in ionizing flux would have produced detectable 
opacity variation in the OVII Ka fine, which instead was not observed. For the z — 0.027 
system, the OVII lines looks stronger in the ACIS-LETG than in the HRC-LETG spectrum 
although their EWs of ~ 2.5 mA and ~ 1.8 mA respectively, are consistent with each others 
within the large ±1 mA la error. Unfortunately for the NVII and NVI lines, the strongest 
lines of this system, a comparison cannot be made, since only the HRC-LETG spectrum 
has enough counts to be sensitive to the detected columns: other HRC-LETG spectra with 
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similar sensitivity limits are needed to strengthen this test. 

We conclude that the WHIM identification is the most likely for both the observed 
z = 0.011 and z = 0.027 systems. 



7.7. Number Density and Mean Cosmological Mass Density of WHIM 

Filaments 

Despite the limited statistics of just two filaments detected along a single line of sight, 
our two detections of WHIM filaments along the fine of sight to Mkn 421 allow us to obtain 
the first crude estimate of the number density of WHIM filaments per unit redshifts above a 
threshold column density, and mean cosmological mass density (Nicastro et al., 2005). The 
expected number density of WHIM filaments above a given ion column density is usually 
computed for the two most abundant ions of the oxygen, OVII and OVIII (e.g. Hellsten 
et al., 1998, Fang, Bryan, & Canizares, 2002). The smallest OVII column density that we 
measure, is Nov// = (7 ± 3) x 10"*^^ cm^^, from the WHIM filament at 2; = 0.027. Using 
this as a threshold, gives a number of OVII WHIM absorbers per unit redshift and with 
Nov// > 7 X 10^^ cm-2, of ZW/A;^ = 2/0.03 ~ 67^^^ (Gehrels, 1986). This estimate of 
dMovii/ is slightly above but consistent with theoretical predictions, which predict ~ 30 
OVII WHIM filaments down to the Nov// > 7 x lO^'' cm'^ threshold and up to z = 0.03 
(Fang, Bryan & Canizares, 2002: from their Fig. 13). Figure 18 shows the predicted curves 
Movii{z) versus z, for three different Nov 11 column density thresholds: 10^^ (solid curve), 
4 X 10^^ (dotted curve) and 7 x 10^^ cm~^ (dashed curve; Fang, Bryan & Canizares, 2002). 
Our data point, Movn{z = 0.03; A'^ov// > 7 x lO"*^^) is also shown with ±lo" errors, and 
is consistent with the predicted curve 

Movn{z]Novn > 7x 10^^): dashed curve in Fie;. 
18. The number density of OVII absorbers with Nov// > 7 x 10^"^ cm~^ is between 1.2 
and 8.3 times that of the corresponding number density of OVI intergalactic absorbers with 
Nov/ ^1.2 x 10^^ cm-2 (i.e. ~ 60 times thinner than the OVII filaments), of dN / dz ^ 19 ±?> 
(Danforth & ShuU, 2005). 

Following Savage et al. (2002) (and approximating the formula to the low-redshift 
regime), we estimate the baryon mass density of the WHIM, in units of the current critical 
density Pc as: 

pc "Mjkn421 

where yU is the average molecular weight (taken to be 1.3), N\j is the equivalent H column 
density of the WHIM filament i and dMknA2i — (zc/Hq) = 128/ifo^ Mpc, is the distance 
to Mkn 421. Averaging on the lower and upper limits of Nh for the two filaments (§7.4, 
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0.01 0.1 1 
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Fig. 18. — Predicted curves Novii{z) versus z, for three different Nqvii column density thresholds: 10^^ 
(solid curve), 4 x 10^^ (dotted curve) and 7 x 10^^ cm~^ (dashed curve; from Fang, Bryan & Canizares, 
2002). The data point, Movii{z = 0.03; Nqvii > 7 x 10^^) is also shown with ±la errors, and is consistent 
with the dashed curve for Afoviiiz; Nqvii > 7 x 10^^). 
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7.5), and including low number statistics errors for a Poissonian distribution (Gehrels, 1986) 
gives fib = 0.027+o;ofg x 10~^^/^^-^. Comparing this value with the fib expected shows that 
virtually all the 'missing' baryons at low redshift are to be found in the WHIM (Nicastro et 
al., 2005). 

8. Summary and Conclusions 

In this paper we presented the first two high significance detections of X-ray WHIM 
filaments, in the very high signal to noise Chandra-LETG spectrum of the blazar Mkn 421 
{z = 0.03). A recently published companion paper (Nicastro et al., 2005) focus on the mea- 
surement of the cosmological mass density in the OVIl WHIM, derived from the detections 
presented here. Our main results are here summarized: 

• We detect weak metal absorption from two intergalactic ionized absorbers along the 
line of sight to Mkn 421, at average redshifts of cz = 3300 ± 300 km s~^ {z = 0.011) 
and cz = 8090 ± 300 km s"^ {z = 0.027). 

• The z = 0.011 X-ray absorber is detected through NVll, OVll and OVlll absorption, 
while the z = 0.027 X-ray absorber is detected through CVI, NVI, NVII, OVII and 
NelX. 

• The conservative significances estimated for the z — 0.011 and z — 0.027 X-ray ab- 
sorbers are 3.5a and 5.8(7. 

• The z — 0.011 system is found at a redshift consistent with a known HI Lya absorber, 
in a sub-dense galaxy environment about 5 times less populated than the average in 
the local Universe. The 3a upper limit temperature of the HI absorber, as derived 
from the observed 3a FWHM of the HI Lya, is inconsistent with the allowed range 
of temperatures of the X-ray absorber. If the two absorbers are related, they imply 
a multiphase IGM. The OVl detection in the FUSE spectrum of Mkn 421 at 2; = 
0.011 ± 0.001 is hampered by a broad detector flaw. 

• The z — 0.027 system lies ~ 13 Mpc from Mkn 421, and is found in an average density 
galaxy environment. No HI Lya or OVI doublet are detected in the GHRS and FUSE 
spectra of Mkn 421 at z = 0.027, down to the 3a sensitivity limits of Nhi > 8.5 x 10^^ 
cm~^ and Nqvi > 1.4 x 10^^ cm~^. 



• We find that a scenario in which the gas is primarily coUisionally ionized, but still 
undergoing additional, and marginal, photoionization by the diffuse UV- X-ray back- 
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ground and the proximity effect of Mkn 421 during bursts, allow us to define ranges of 
temperatures and metallicity ratios for both intervening systems. 

• For the z — 0.011 absorber we find a temperature range T — (0.6—2.5) x 10^ K. A rather 
tight lower limit can be set on the metaUicity of this IGM system, of [0/H]> —1.47. 
Only a weak lower hmit can be set on the relative N/0 ratio, of [N/0]> —2.05, 
consistent with the Solar value. Given the allowed range of temperature for the z — 
0.011 absorber, equivalent H column densities of N/^ = [(1.7 ± 0.7)T6 — (0.1 ± 0.6)] x 
10^^ X 10^['^/^l-i (with [N/0] = 0) are implied. These in turn imply a thickness of the 
absorber along our line of sight oiD = [(0.55±0.23)T6- (0.03±0.19)] x \^-^^ I ^^-\n^)-_\ 
Mpc. 

• For the z — 0.027 absorber, we find that solution exist in the hybrid equilibrium 
collisional/photoionization scenario and for = 10^^ cm^"^, only for a narrow (AT ~ 
2 X 10^ K) range of temperature around T= 1.4 x 10^ K, and only for rather high N/0 
ratios compared to Solar: [N/0]> 0.65, and moderately high metallicity [0/H]> —1.32. 
The breadth of the AT interval, around the T= 1.4 x 10^ K solution, depends on the 
metallicity ratio: values of [0/H] larger than -1 will broaden the allowed range of 
temperatures. For [N/0] = 0.65 we find Njj = [(1.67 ±0.76)T6 + (0.71^1.15)] x lO^'^ x 
10-[o/i/]-i cni-2. This implies D = [(0.53 ± 0.24)T6 + (0.23 0.37)] x 10-[°/^l-i(n6)l^ 
Mpc, assuming constant baryon density. 

• All the physical properties of these two intervening absorbers are consistent with theo- 
retical predictions for the WHIM. We therefore identify these two intervening absorbers 
with two tenuous WHIM filaments absorbing the X-ray radiation of Mkn 421 along 
our hue of sight. This represents the first > 3(7 detection of WHIM filaments in the 
X-rays. 

• Based on our detections we set for the first time an estimate of the number density 
of OVII WHIM filaments with ^ovii > 7 x lO^'' cm^^ in the local Universe, of 
dM /dz — 67;tl3, (1-2-8.3) times that of the corresponding, but ~ 60 times thinner, 
OVI absorbers. Finally, from our equivalent hydrogen column density measurements, 
we derive a mean cosmological mass in X-ray WHIM filaments, in units of the current 
critical density of the Universe, of Vli, = 0.027+0 019 x 10'^^/^^-^ consistent with both 
model predictions and the estimated number of 'missing' baryons at low redshift. 

This work only begins the study of the Warm-Hot IGM outside the Local group. The 
large number of ion species showing up in the X-ray spectra imply a rich field of investigation. 
Higher significance detections, along several other lines of sight, and at higher redshift, of 
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He-like and H-like transitions from C, N, O and Ne, will allow us to measure the relative 
metallicity ratios at different redshifts. This will fundamentally contribute to the assessment 
of the still very poor measurements of metal production at relatively high redshift, and will 
help refining models for galaxy-IGM feedback, possibly providing invaluable hints for the 
main mechanisms responsible for this feedback (i.e. galaxy superwinds, versus quasars winds, 
etc.). Moreover, the potential for the simultaneous detection, in the same band, of inner-shell 
transition from lower ionized species of the same ions, will allow us to distinguish between 
different ionization scenarios, and to assess the importance and frequency of multiphase IGM. 
Reconstructing the distribution functions of several ions in the WHIM (i.e. dj\f /dzdN), will 
also probe very useful in distinguishing different cosmologies (Fang &: Canizares, 2000), and 
possibly tracing the distribution of dark-matter potential wells in the local Universe. 

This work clearly shows that an extremely efficient way to detecting the WHIM with 
Chandra and A^eu;ton-XMM, is observing sources while in bright outburst phases. Mkn 421 
is unique: no other extragalactic sources in the sky (except Gamma-Ray Burst X-ray after- 
glows, e.g. Fiore et al., 2000) reach the level of Mkn 421 during the observations presented 
here. However, as shown in Fig. 18, the expected number of OVII WHIM filament down 
to a given Noy// threshold, increases linearly with redshift. About 2 OVII filaments with 
Noy/7 > 10^^ cm~^ are expected up to 2; = 0.5 along a random line of sight. Such OVII 
columns can easily be detected with 100 ks Chandra and/or Newton-XMM. observations of 
background sources flaring at a level even 10-50 times lower than that reached by Mkn 421 
during the observations presented here, or 100-500 ks observations of pre-selected 2; ^0.3 
blazars with quiescent soft X-ray flux level of Fo.5_2feey ~ 1 — 2 mCrab. 

A much greater potential is offered by the large collecting area of the spectrometers 
on the two planned X-ray missions, Constellation-X and XEUS . Constellation-X will 
allow us to detect WHIM systems as weak as those presented in this work, for possibly 
tens of different lines of sight, and thicker filaments will be detected towards hundreds of 
different sightlines. A resolving power of R~ 500 at 20 A will not allow us, m many 
cases, to resolve the lines and detect multiphase structures in the IGM. This will in turn 
hamper our capabilities of constraining the physical properties of the baryons in this diffuse 
web of WHIM and so of measuring fl^^^^ . A resolving power R^SOOO — 5000 is needed 
to resolve thermally broadened OVII fines in gas at T = 10^ K and so to (a) disentangle 
the Doppler parameter b from the ion column density, and so derive the saturation-corrected 

^^http: / / constellation.gsfc.nasa.gov/ science/index. html 
-•^^http: / /www.rssd.esa.int/index.php?project=XEUS 

•^"http: / / constellation.gsfc.nasa.gov/ science/design/spectraljresolution.html 
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contribution to Qb (b) measure the intrinsic internal turbulence in the filaments, and compare 
it to the galaxy velocity in the field, (c) allow for separation of multi-phase components, and 
in general, for a clear physical diagnostics (Elvis & Fiore & the Pharos Team, 2003). 
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